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Abstract 

d ' We examine the problem of two particles confined in an isotropic harmonic trap, which interact 

via a finite-ranged Gaussian-shaped potential in two spatial dimensions. We derive an approxi- 

^: 

^3 ! mative transcendental equation for the energy and study the resulting spectrum as a function of 

^-^ ■ the interparticle interaction strength. Both the attractive and repulsive systems are analyzed. We 

study the impact of the potential's range on the ground-state energy. Complementary, we also 
explicitly verify by a variational treatment that in the zero-range limit the positive delta potential 
in two dimensions only reproduces the non-interacting results, if the Hilbert space in not truncated. 
Finally, we establish and discuss the connection between our finite-range treatment and regularized 
zero-range results from the literature. 

PACS numbers: 03.75.Hh, 34.20.-b, 03.65.-w 



1 



I. INTRODUCTION 



In recent years there has been an increasing interest in two dimensional (2D) quantum 
systems. The condensed matter community has long been investigating 2D quantum effects, 



for example, in relation to superfiuid films 



fla, 



high-temperature superconductivity 



3, arid low-dimensional materials, sueh as graphene 0, Q. The rapid progress in atomic 
trapping and cooling now also allows to study quantum systems with reduced dimensionality 
in the context of ultracold trapped gases [8]. Degenerate quasi-2D Bose and Fermi gases 
have already been produced in highly anisotropic "pancake" traps 9l-|l4|. This opens up 
the unique possibility to investigate the rich palette of 2D quantum effects and phases in a 
highly controlled environment. 

In dilute systems, such as the ultracold trapped gases, interactions are usually described 
by a two-particle zero-range effective potential. This approach has been particularly fruitful 
in one dimension, where the delta function interaction is well behaved and there is a simple 
relation between the scattering length and the interaction parameter. In two and three 



15 



17 



dimensions, however, the delta function is not a self-adjoint operator [15| 
rise to various anomalies 



16|, which gives 



21] . In particular, there is no scattering from a positive 2D 



delta function and, for a negative one, the bound-state energy diverges. Several approaches 
to overcome the arising problems and to model zero-range interactions in 2D have been 
proposed in the literature: self-adjoint extensions 22,|23[, renormalization techniques 



9 

261 



regularization of the delta potential 



16 



24 



25|, and modified boundary conditions 



The analytically tractable problem of two harmonically trapped particles interacting via 
a zero-range potential in 2D has previously been addressed using a regularized delta function 
25 | and Bethe-Peierls boundary conditions j^^). These works have obtained and studied 



16 



the spectrum of the particles with zero-range interaction as a function of the 2D scattering 
length 



16 



In the present work we examine the problem of two harmonically trapped particles in 
2D interacting via a finite-ranged two-body potential modeled by a Gaussian function. We 
derive an approximative, yet accurate, transcendental equation for the energy, and present 
the resulting spectrum. In particular, we study the energy levels for both positive and 
negative interaction and, furthermore, explore the role of the range of the potential on the 
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ground-state energy. Finally, we establish and discuss the connection between our finite- 
range results and previous zero-range works. 

The article is organized as follows: in section II we derive the general secular equation 
for the energy of two trapped particles interacting via Gaussian-shaped two-body potential 
in 2D. In section III, by utilizing a variational treatment, we briefly discuss the limit of 
a positive non-regularized delta potential and show that the energy spectrum of the non- 
interacting system is altered only as a consequence of truncating the Hilbert space. In 
Section IV A we present an efficient high-performance approximation for the finite-range 
interaction, and derive an equation for the energy of the two particles. Then, in section 
IV B, we study the resulting energy spectrum and in section V compare our finite-range 
findings to zero-range results from the literature. Finally, in section VI, we summarize our 
results. Supplemental derivations and numerics are deterred to the appendices. 



II. THE EIGENVALUE EQUATION 

We consider two particles in an isotropic harmonic trap, which are interacting via a 
normalized two-body Gaussian-shaped potential. The Hamiltonian of the system is 



i=l ^ ' 



(1) 



V{fi-f2) = —^e (2) 

Here V is the 2D Nabla operator and f = (x, y). The problem can be separated into a non- 
interacting center of mass and an interacting relative part. With the standard definitions 
of a reduced mass, \x — m/2, and a total mass, M — 2m, as well as center of mass and 
relative coordinates R = \{fi + f2) and f = fi — f2, the Hamiltonian can be rewritten as 
H = Hem + Hrei with: 

^ + ^Ma;2i?^ (3) 

Hrei = -^V^ + ^i^uj\^ + XoV{r). (4) 
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Here Hem is the Hamiltonian of the 2D quantum harmonic oscillator whose solutions are 
well known. From now on we concentrate on the relative part which we further write as 
Hrei = Hq + \oV{r), where Hq reads 

H. ^ -fvi + (5) 

We start by constructing a solution of the time-independent Schrodinger equation 
from the eigenstates of Hq. We take only eigenstates with zero angular momentum, which 
we denote ipk{f^)- The respective energies are Ek = {2k + 1)^. We note that all states with 
non-zero angular momentum vanish at the origin and, therefore, should not be significantly 
perturbed by the Gaussian potential if its width is sufficiently small. After substituting the 
expansion \1/ = Xli^o '^'^V'A: into the Schrodinger equation H^ei^ = , and after projecting 
onto a state (fki{r), we arrive at the secular equation 

Ck'iEk' - ^) + Ao ^ Cfc / ipl,{r)V{r)^k{r)df = 0. (6) 

The integration is taken over the whole 2D plane with df = 2'Krdr. To proceed further we 
need to evaluate the matrix elements appearing in the above sum. As we show in Appendix 
A, this can be done analytically. The result of the calculation is 

/ X k' + k + l 



Here, 2F1 is the Gauss hypergeometric function |28| and / = y ^ is the harmonic oscillator 
length. 



III. SOLUTION FOR CONTACT POTENTIAL 



Before proceeding to the results for a finite-ranged interaction, let us first examine the 
limit s — 7- 0, in which case the normalized Gaussian-shaped potential defined in Eq.([2]) goes 
into a delta function. For s = Eq.([7]) takes on the form Ik',kiO) = and the matrix 
elements are independent of the indices k' and k. In this limit Eq.(l6]) reduces to 

00 ^ 

A:=0 
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The sum appearing in the above equation runs over all indices k. After rearranging the 
expansion coefficients we obtain 



Ck' = ^, (9) 

where C is a parameter which can depend on the energy E, but is the same for all k'. By 
substituting the above expression for the coefficients back into Eq.® and dividing by 
nl^/Xo, we arrive at 

— + g 2t + l-g/fi^ =°' (1°) 



29|. The 



The sum in the above equation is a general harmonic series and is divergent 
anomaly associated with a contact potential in 2D stems from this diverging sum in the 
current treatment. To obtain a meaningful expression we ffist truncate the sum at finite 
and then examine the behavior for N ^ oo 

hojnP ^ 1 _ 

^ ^2k + l- E/hw ~ ^ ^ 

k=0 

From the above equation we can determine the energy spectrum of the two trapped particles 
interacting via a delta potential in 2D for a given truncation A^. 

Let us consider the solution of Eq. ( fTTl) closest to one of the poles appearing in the equa- 
tion, say, the pole specified by A; = k'. We may write E/hw = 2k'+l+Ak'. This immediately 
leads to 

1 hwrrl'^ ^ 



Ak' Ao ^ 2k + l- E/hw' 

kj^k' 

For large A^ the sum on the right-hand side (rhs) grows logarithmically with A^ for any 
value of E ^ 2k + 1 {k ^ k'), and hence Ak' approaches zero as ~ 1/lnA^ . Thus, by 
increasing A^ we can make the energy levels arbitrarily close to the eigenvalues of the respec- 
tive non-interacting system. These considerations explicitly show that in 2D the positive 
(non-regularized) delta potential modifies the spectrum of two trapped particles only as a 
consequence of restricting the Hilbert space. We point out that the non-interacting values 
are approached logarithmically because the series in Eq. lfTT]) diverges logarithmically. We 
illustrate this in Fig.([T]), where we plot the energy of the lowest state obtained from Eq. (fTT|) 
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for increasing A^. We stress that the above considerations are rigorous and Eq. lfTT]) can be 
viewed cis db variational ansatz (see Appendix B). 




1 I I I I I I I 

1 10^ 10^ 10^ 10^^ 10^^ 10^^ 

N 

Figure 1: Ground-state energy (in the center of mass frame) of two particles in a harmonic trap 
interacting via a (non-regularized) delta potential in 2D as a function of the size of the Hilbert 
space. The results are obtained by numerically solving Eq. fjlip with Aq = 1, / = 1 and different A'^. 
Notice the logarithmic scale on the x-axis. All quantities are dimensionless. 



IV. FINITE-RANGE POTENTIAL 

A. Efficient high-performance approximation 

In order to obtain an (approximative) equation for the energy of the two particles in the 
finite-range case, s > 0, we proceed analogously to the treatment above and make an ansatz 
for the expansion coefficients Ck'- To this end we use an ansatz similar, but not identical, to 



6 



the one obtained in the Brillouin-Wigner perturbation theory (see discussion in Appendix 
C) 



c. = zMm^c (13) 

and obtain an equation for E by substituting this expression into Eq.(l6]). Since we have 
-^o,fc(0) = = const., this ansatz ensures that for s — > we exactly recover the delta 
potential limit of Eq.([9]), which was discussed in the previous section. We stress that the 
current approximation in not variational, in contrast to the case s = 0. However, we 
establish the high accuracy of the treatment by comparing to full direct diagonalization and 
zero-range results (see Appendix D and section V). 

By substituting (1131) into Eq.([2]) and dividing by Aq, we obtain 

J^'^'^' + ^ Sk/fu. - E/fu. = °- 
For each index value k' the above expression gives an equation for E. By setting fc' = 
the matrix elements Ik',k{s) take on the form Jo,fc = i J equation for E 

becomes 

-^-(0,0 + V 77 ETT^ = 0- 

^0 f^^Ek/hw - E/fuj 

P 

The series y^^g ^^i^^'^^iy^ can be expressed in terms of the Lerch transcendent function 
$(;2,s,a) IsO] 



^ Eklhw - Ejhuj (7r/2)2 2A; + 1 - E/hw \ + 2^^^ (1 + {s/iyf 

(16) 

The final equation for the energy of the two trapped particles in the center of mass frame 
reads 

^ V(i+(./o^f' 2 J _hw 

27r/2 (1 + (s//)2) Ao' ^ ^ 
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Eq. (fT7|) is the main analytical result of the paper. It allows one to obtain the energy 
spectrum of two trapped particles interacting via a Gaussian-shaped potential in 2D for 
given parameters Aq, s and /. 



B. Energy spectrum 

In this section we present the energy spectrum (in the center of mass frame) of two 
trapped particles interacting via a Gaussian-shaped potential in 2D, which results from 
Eq. (fT7|) . From now on we fix the harmonic oscillator length to Z = 1. We first set s/l = 0.1 
and explore the dependence of the energy levels on the parameter Aq. The energies of the 
first three states versus Aq are plotted in Fig. ([2]). For Aq > the energies are always above 
the respective non-interacting values and the system is repulsive. As expected, at Aq — )■ 
we recover the non-interacting values E/hu = 2k + 1, which correspond to the poles of the 
left-hand side (Ihs) of Eq.( fTT|) . For Aq < the particles are interacting attractively and the 
ground state energy quickly diverges to ~oo for Aq — > — oo. 

The Gaussian-shaped two-body potential also allows us to study the role of the range, s, of 
the interaction. In Fig. ([3]) we show the dependence of the ground-state energy on s for three 
different values of the parameter Aq. For repulsive interactions, i.e., Aq > 0, we find that 
with decreasing s the ground-state energy logarithmically approaches the non-interacting 
value huj. This is in agreement with our discussion of the 2D delta potential in section UTTl 
We note that this behavior is also in agreement with the formal result in 3ll, where it is 



proven that in two and more dimensions the solutions of the Schrodinger equation are not 
affected by positive potentials with vanishing support. For attractive interactions (Aq < 0) 
the dependence on s is more pronounced and we observe that in the limit s — ?■ the energy 
of the ground state diverges to — oo for any negative Aq. This result is also consistent with 
previous studies, where it was observed that the attractive (non-regularized) delta potential 
yields a bound state with an infinitely negative energy. This is, in fact, the starting point 



for renormalization treatments 



19 



32|. 
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Figure 2: The energies of the first three states (in the center of mass frame) of two trapped particles 
in 2D interacting via a Gaussian-shaped potential versus the interaction strength Aq- The dashed- 
dotted lines show the energies of the respective non-interacting system. Here / = 1 and s/l = 0.1. 
All quantities are dimensionless. 

V. COMPARISON WITH ZERO-RANGE RESULTS 

The energy spectrum of two harmonically trapped particles with a zero-range interaction 



in 2D has previously been obtained in the literature using a regularized delta potential 



16|, 



general scattering arguments 25|, and modified boundary conditions [27]. In all three works 
the authors derive a transcendental equation for the energy, which can be written in the 
form 

Here 'il){x) denotes the Digamma function, a2D is the 2D scattering length, and I is the 
harmonic oscillator length. The interaction strength is controlled through an interaction 
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Figure 3: Energy of the ground state (in the center of mass frame) of two trapped particles in 
2D interacting via a Gaussian-shaped potential versus the width of the interaction, s, for different 
choices of Xq/Ttw. See legend inside the graph. Notice the logarithmic scale on the x-axis. All 
quantities are dimensionless. 



parameter defined as 



within the works 



16 



n 



25 



16 



25] . The constant factor A, which is slightly different 



27| . depends on the exact form for t 



le effective range expansion 



25| ) and has no consequence 



that the respective authors use (see the discussion at the end of 
on the following analysis. 

In order to relate our finite-range result for the energy of two trapped interacting particles 
to the zero-range treatments, we employ a Taylor expansion of the Ihs of Eq. (fT7|) around 







33| 



1 l-E/hw 





(i+(wr' 



27r/2 



i-E/njuj\ 

2 J 



+ ln (^)+7 



+ 0[s]. (19) 



Here 7 ~ 0.577(2) is the Euler-Mascheroni constant and ip{x) is again the Digamma function. 
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By neglecting terms of order 0[s] and after rearranging, we can rewrite Eq. (fT7|) for s as 



~ fl-E/nw\ 277/2 

with A{s) = — In ~7- Evidently, the above equation for the energy has the same form 

as the literature result in ( fTSl) . and by an appropriate choice of Aq and s our ansatz reproduces 
the zero-range spectrum. By comparing the rhs of Eq. (120]) and f llSp we immediately see that 
for weak interactions) the interaction parameter is equal to in our finite-range analysis 
34j |. Thus, the interaction parameter in 2D is proportional to the factor Aq in front of the 
two-body potential [see Eq.([T])]. 

The Taylor expansion in ( IT9|) leads to two significant differences between the zero-range 
spectrum and the finite-range spectrum already presented in Fig. ([2]). First, when the inter- 
action is repulsive, i.e., Aq > 0, Eq. f|20p yields an additional deeply-bound state which is not 

35] for the case of hard spheres in 



present in the original finite-range spectrum (see also 
three dimensions, where also a redundant state appears in the zero-range pseudopotential 
approximation). This state appears due to the different asymptotic behavior of both sides 
of (fT9|) in the limit E/fm — oo. While the Ihs of (fT9|) converges to 0, the expanded rhs 
diverges to -|-oo, see Fig.(j4]). This also leads to a second difference. The finite-range result 
yields a bound state with an energy approaching — oo for Aq — )■ — oo, while the zero-range 
equation yields a finite value which corresponds to the zero crossing appearing for negative 
E/hw in Fig. ([2]). This demonstrates that the short-ranged 2D Gaussian-shaped interaction 
potential and the zero-range potential approximate each other well except for the ground 
state. 

Equating the rhs of Eq. (l20l) and (fTSl) allows us to obtain a connection between the 2D 
scattering length, a2D, and the parameters of the Gaussian-shaped potential Aq and s. Using 
A = 2 ln(2) — 27 from j25|, and noticing that the authors use the definition of the harmonic 
oscillator length with the full mass, we obtain 

a2D ~ v^se"^-A^ (21) 

for the 2D scattering length of a Gaussian-shaped potential of width s/l ^ 1. We find 
good quantitative agreement when we compare the above expression for a2D with numerical 
values from [36]. For a = 2~^^'^s = 0.1 and g = Xq = 1 the relative difference between our 
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full 

expansion 




Figure 4: Plot of the Ihs (thick hne) and the rhs (dashed hne) of Eq. (|19p versus Ejhio for / = 1 
and sjl = 0.1. For E/fiu] > — 1 the two curves are practically indistinguishable. The vertical 
dashed-dotted lines show the poles of both sides of Eq. (|19p . Notice the logarithmic scale on the x 
axis. See text for discussion. All quantities are dimensionless. 



analytical estimate and the numerical value is about one percent. Even in the regime of 
a = = 1 and g = Xq = 10, where the accuracy of our treatment is limited, the relative 

difference in only about 8 percent. 



VI. SUMMARY 

We analyzed the problem of two harmonically trapped particles in 2D, which interact 
via a finite-ranged Gaussian-shaped two-body potential. We derived an approximative tran- 
scendental equation for the energy which is demonstrated to be highly accurate for the 
ground state. Moreover, our ansatz works well also for the excited states (see Appendix 
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D). Using the Gaussian-shaped potential we were able to directly study the dependence of 
the ground-state energy on the range of the interaction. We found that the effect of the 
short-ranged interaction on the ground-state energy vanishes logarithmically with deceasing 
of the potential range, s, for all positive interaction strengths. This study is complemented 
by a variational treatment which shows that in the limit of a (non-regularized) delta poten- 
tial, i.e. s — )■ 0, the energy spectrum of the two (repelling) particles can be made arbitrarily 
close to the respective non-interacting one by increasing the size of the Hilbert space. Fi- 
nally, we established and discussed the connection between our finite-range result and earlier 
zero-range treatments reported in the literature. 

Our analysis shows that a Gaussian-shaped two-body potential and its zero-range pseu- 
dopotential give similar results for both repulsive and attractive interactions except for the 
lowest eigenstate of the latter. Here, for repulsion the zero-range pseudopotential leads to an 
additional dimer bound state which is not connected to the unperturbed system when the 
interaction is switched off. Consequently, this dimer state is difficult to reach when the in- 
teraction is smoothly switched on. When going beyond two particles, analytical treatments 
quickly become inaccessible, leaving one only with a numerical recourse. For the numerical 
many-body simulations one usually prefers to use smooth, finite-ranged, model interaction 
potentials. Our results support the applicability of the Gaussian-shaped potential as a model 
of the two-body interaction in such simulations in two dimensions, see for example [36]. 
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Appendix A: Evaluation of the matrix elements 

Here we outline the calculation of the matrix elements in Eq.(j7]) of the main text. The 
radial part of the zero angular momentum eigenfunctions of the 2D harmonic oscillator can 
)e written in terms of the Laguerre polynomials Lk{x). With ^ = y ~ the eigenstates read 



37| 



„2 



Mr) = V2l-'e-^Lk ( ) . (Al) 



We need to evaluate the integral 



-r2 







Ik',k = 2/-' / e-^Lk, C-^] ^ e-^Lk (^]rdr. (A2) 



We start by setting A = 1 + 4- and rewriting 



= ^ r e^^^^fc' 1 L,(^-^)r dr. (A3) 



TTS^ 







P "\P 



Now we define a new variable p = jjr"^ and obtain 
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1 /""^ 

h'^k = -r-^ / e-^^Lk> (A-V) Lj, (A-^) dp. (A4) 
To get rid of the A~^ factor in the argument of L{A~^p) we use the known multiphcation 



formula for Laguerre polynomials 



28 



L,(Ax) = 5^(^)a^(1-A)'=-L, (x). 



(A5) 



By substituting the above expression for Lk{A~^p) and Lkr{A~^p) into the integral of 
Eq. (lA4p . and by making use of the orthogonality of the Laguerre polynomials with respect 
to the weight function e~P, we obtain 

min{k,k'} , . , . 

^'■.-^ E (*)f:)A-(l-A-')--. (A6) 

We can write the above expression as 

In order to express the sum in terms of the Gauss hypergeometric function 2Fi(a, 6; c, z), we 



rewrite the above expression using the Pochhammer symbols {x)i for the rising factorial 



28 



/ . -. \ k+k' minfk.k'} , , , , , ,n ( 1 

^'^''^-^l^j ^ — ^ — • 



1=0 



We can now directly use the definition of 2Fi(a, 6; c, |28| and obtain 



1 ( „ / ., . , / 1 ^ ^ 



Finally, by changing back to the original variables, we arrive at Eq.([7]) of the main text. 



Appendix B: Variational energy 

Here we show that the approach in Section II A for the energy spectrum of the non- 
regularized delta potential is essentially a variational treatment with = Xlilo '^fcV^fc and 
Cfc = ■ For this purpose, let us consider E just as a parameter which fulfills Eq. lfTT]) for 
a given Aq and fixed finite N . The variational energy then reads 
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\k=0 / \k'=0 

\ ^ 1 

const. 



N 

E 

fc=0 



const. 
Sk-E 



Ek — E ' -nP ^ 6k' — E 

. k'=0 



N 



k=0 



const. 
Ek — E 



^k 



E. 



We thus see that the variational energy is indeed equal \,o E. In the last step we used the 
normalization condition for the wavef unction XlfcLo kfcp = 1- 



Appendix C: Ansatz for Ck in the finite-range case 

Here we discuss the ansatz for the expansion coefficients in Eq. (fT^ . We stress that the 
approach is general. Let us consider a Hamilton operator 

H = Ho + XoW (CI) 

and assume that the eigenstates and eigenenergies of Hq are known. We denote them by 
with ifo|0fc) = £k\4>k)- For an eigenstate |\1') = '^Ck\(f)k) of the full Hamilton operator 
H we can write 



{Ho + XoW)\^) = E\^), 



or equivalently 



(E - Ho)\<iJ) = \oW\<iJ). (C2) 
After projection on an unperturbed state (0^1 we express the above equation in the form 



{E-ek){<pk\'^) = Xo{<pk\W\^). 
From here we obtain a symbolic expression for the expansion coefficients 



(C3) 



Ck 



Xo{(pk\W\^) Xo{(pk\W\^) 



(E-Sk) 



{sk - E) 



(C4) 
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Until now we have not used any approximations. Of course, the above expression is only 
an implicit one, because the unknown state |\E'), which is itself dependent on Cfe, appears 
on the rhs. Our ansatz in Eq. (fT3|) of the main text consists in taking only a "first-order 
approximation" for the expansion coefficients, i.e. we substitute |^) = {(po) in (1C4I) and 
obtain 

= - (e,-E) ■ e^^' 
Of course, we can also substitute |^) = l^^')' k' ^ into Eq. OCSP and obtain similar 
expressions which might be more efficient when studying excited states. 



Appendix D: Numerical comparison 

In this section we compare energies obtained by solving Eq. ( fT7|) with the numerical values 
determined by the full direct diagonalization of the Hamiltonian using the matrix elements 
in Eq.(I7]). In Tab.(ITl) we show the relative difference between both values, i.e. , 
where denotes the value obtained from direct diagonalization and E denotes the energy 
obtained by solving Eq. (fT7|) . The comparison is performed for the ground state (GS) and 
several excited states with zero angular momentum [see the first column in Tab.([T])]. We 
observe an excellent agreement between the values from both methods, with in the 

order of 10 ~^ and below, see Tab.(Il]). 

We also perform a study for a fixed width, s = 0.2, and variable Aq. The results are shown 
in Fig. ([5]), where we plot the values of the ground-state energy as obtained from Eq. (fT7|) and 
direct diagonalization. The values obtained from Eq.f ll7p are in excellent agreement with 
the ones from direct diagonalization for Xo/hw up to 4. Analogous behavior was also found 
for the excited states. 
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0.5 1 1.5 2 2.5 3 3.5 4 

^0 / 



Figure 5: Comparison between the ground-state energy £"0 versus Aq obtained from Eq. fjlZp (thick 
hue) and that computed by direct diagonahzation (crosses). The width is fixed to s = 0.2 and 
I = 1. See text for more details. All quantities are dimensionless 





s = 0.1 


s = 0.2 


s = 0.4 


s = 0.5 


GS 


6.1 X 10""^ 


7.4 X 10"^ 


6.9 X 10"^ 


5.8 X 10-^ 


1st 


2.1 X 10""^ 


9.2 X 10"^ 


4.9 X 10"^ 


1.7 X 10-3 


2nd 


1.1 X 10"^ 


1.6 X 10"^ 


1.8 X 10-3 


4.1 X 10-3 


4th 


3.4 X 10"^ 


1.5 X 10"^ 


3.5 X 10-3 


5.2 X 10-3 


8th 


4.7 X 10-6 


5.7 X 10-"^ 


3.3 X 10-3 


3.4 X 10-3 


16th 


3.9 X 10"^ 


1.0 X 10-3 


1.7 X 10-3 


1.4 X 10-3 



Table I: Numerical values of for the ground state and several excited states with Aq = 1.0 

and four different choices of s and 1=1. E^i^i is the value obtained by direct diagonahzation of 
the Hamiltonian and E is computed using Eq. (|17p . 
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